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ON . ABSTRACT 

o ■ 

IL/ ' Aims. Galaxy clusters are key targets in the search for ultra high energy particle accelerators. The Coma cluster represents one of the best candidates 
. 5^ ' for such a search owing to its high mass, proximity, and the established non-thermal radio emission centred on the cluster core. 

Methods. The H.E.S.S. (High Energy Stereoscopic System) telescopes observed Coma for ~8 hr in a search for y-ray emission at energies > 1 TeV. 
j_j The large 3.5° FWHM field of view of H.E.S.S. is ideal for viewing a range of targets at various sizes including the Coma cluster core, the radio- 
ed relic (1253+275) and merger/infall (NGC 4839) regions to the southwest, and features greater than 1° away. 

Results. No evidence for point-like nor extended TeV y-ray emission was found and upper limits to the TeV flux F(E) for E > 1, > 5, and 

> 10 TeV were set for the Coma core and other regions. Converting these limits to an energy flux E^F(E) the lowest or most constraining is the 

E > 5 TeV upper limit for the Coma core (0.2° radius) at ~8% Crab flux units or ~ 10^'^ph cm"- s"'. 

Conclusions. The upper limits for the Coma core were compared with a prediction for the y-ray emission from proton-proton interactions, the level 
of which ultimately scales with the mass of the Coma cluster A direct constraint using our most stringent limit for E >5 TeV, on the total energy 
content in non-thermal protons with injection energy spectrum oc E^~' and spatial distribution following the thermal gas in the cluster, is found to 
be ~0.2 times the thermal energy, or ~ 10''^ erg. The E >5 TeV y-ray threshold in this case corresponds to cosmic-ray proton energies S50 TeV. 
Our upper limits rule out the most optimistic theoretical models for gamma ray emission from clusters and complement radio observations which 
constrain the cosmic ray content in clusters at significantly lower proton energies, subject to assumptions on the magnetic field strength. 
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1. Introduction 

Clusters of galaxies represent the largest gravitationally bound 
objects in the Universe and are thought to be ideal sites for 
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the acceleration of particles. The very long confinement time 
(of order the Hubble time) of the accelerated particles (see 
e.g. Yolk et al. 119961 Berezinsky et al. |1997| l would allow in- 
teractions of the particles with ambient matter and radiation 
fields to produce non-thermal emission from radio to TeV y- 
ray energies. Particles are thought to be accelerated at large- 
scale shocks associated with accretion and merger processes 
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(see e.g. Colafrancesco et al. 119981 Ryu et al. I2003I I. in super- 
nova remnants and galactic-scale winds (Volk et al. I1996I I. tur- 
bulent re-acceleration (Brunetti & Blasi I2005I I and dark matter 
annihilation (e.g. Colafrancesco et al. 120061) . In addition, parti- 
cles may be re-distributed/injected throughout the cluster vol- 
ume via AGN cluster members (EnBlin et al. 1997 Aharonian 
12002! Hinton et al. 120071) . The non-thermal radio emission ob- 
served in recent years from several galaxy clusters (Giovannini 
et al. 119931 Feretti et al. I2004I I represents clear evidence for rel- 
ativistic particle populations in such objects. Further evidence is 
provided by possible non-thermal X-rays observed from a few 
clusters (Rephaeli & Gruber 12002! Fusco-Femiano et al. !2004! 
Eckertetal. "20021). 

Gamma-ray emission in galaxy clusters may come from sev- 
eral processes (see review by Blasi et al. !2007l ). The collision of 
relativistic cosmic -ray (CR) protons with thermal nuclei com- 
prising the intra-cluster medium (ICM) may lead to y-ray emis- 
sion via the decay of neutral pions (Dennison !1980l l. In this con- 
text the fraction 77 of thermal energy in the cluster volume in the 
form of relativistic non-thermal particles is an important parame- 
ter that can determine the level of y-ray emission expected. Since 
the thermal energy content is a function of the cluster mass, 
the most massive and nearby clusters present the best oppor- 
tunity to probe for such y-ray emission. Ultra-relativistic elec- 
trons can also up-scatter target photons such as the cosmic mi- 
crowave background (CMB), infrared, starlight, and other soft 
photon fields) to TeV y-ray energies (Atoyan & Volk !2000! 
Gabici et al. 2003] 120041 ). Given that galaxy clusters may ac- 
celerate particles to ultra high energies (UHE) > lO'^ eV (e.g. 
Hillas !1984! Kang et al. 119961) . y-ray production from inverse- 
Compton scattering by secondary electrons generated when a 
UHE proton interacts with a CMB photon in the Bethe-Heitler 
process (pj — > e^e^ + p') may also result (Inoue et al. 12005! 
Kelner & Aharonian !2008! ). Dark matter annihilation has also 
been considered as a y-ray production channel (e.g. neutralino 
annihilation by Colafrancesco et al. !2006!) . 

Earlier observations in the MeV to GeV y-ray band with 
EGRET only found upper limits only for several clusters 
(Reimer et al. '20031) including the Coma cluster. At TeV ener- 
gies, upper limits (Perkins et al. !2006l ) have been reported for the 
Perseus and Abell 2029 clusters using the single-dish Whipple 
telescope. The most recent TeV observations with stereoscopic 
instruments suc h as H.E.S.S. (Abell 496 and Abell 85 — 
Aharonian et al. l2008b . and with VERITAS (Coma — Perkins 
et al. !2008!) . revealed also upper limits. This work focuses on 
H.E.S.S. observations of the Coma cluster 

Coma (ACO 1656) is one of the nearest (z=0.023) and 
best-studied galaxy clusters. Extended (several arcminutes in 
scale) hard X-ray emission (with so far weak evidence for 
a non-thermal component) has been observed (Rephaeli & 
Gruber 2002, Fusco-Feminano et al. 2004 , Rossetti et al. l2004l 
Lutovinov et al. .2008. Ajello et al. ,2009) , as well as a promi- 
nent non-thermal radio halo (Giovannini et al. !1993! Thierbach 
et al. !2003! ). The latter is clear evidence for particle acceleration. 
Being one of the most massive (M ~ lO'^M©) and nearby clus- 
ters, with detailed multiwavelength observations ranging from 
low frequency radio wavelength to y-rays, the Coma cluster has 
always been considered as the prototypical cluster also for very 
high energy y-ray studies. The Coma cluster is located in the 
northern hemisphere and is visible by H.E.S.S. at moderately 
high zenith angles (average value ~ 50°), which leads to a rela- 
tively high energy threshold (defined as the peak detection rate 
for an ' power-law spectrum of y-rays) of > 1 TeV. Since 
the y-ray spectrum from clusters is expected to be hard and ex- 



tend beyond 10 TeV (basically limited only by the absorption of 
y-ray photons in the cosmic infrared background. At the Coma 
cluster distance, an optical depth of unity is reached for energies 
E ~ 10 to 20 TeV), the energy threshold does not constitute a 
serious problem for our investigation. 

2. H.E.S.S. observations and analysis 

Operating in the Southern Hemisphere, H.E.S.S. consists of four 
identical 13 m diameter Cherenkov telescopes (Bernlohr et al. 
|2003 ). H.E.S.S. employs the stereoscopic imaging atmospheric 
Cherenkov technique, and is sensitive to y-rays above an en- 
ergy threshold of ~0. 1 TeV (Hinton et al. 12004] ) for observa- 
tions at zenith. An angular resolution of 5 to 6' on an event-by- 
event basis is achieved, and the large field of view (FoV) with 
FWHM ~ 3.5° (Aharonian et al. '2006b) permits survey cover- 
age in a single pointing. A point source sensitivity of ~1% Crab 
flux (~ lO"'-' erg cm^^s"' at 1 TeV) is achieved for a 5cr 
detection after ~25 hr observation. Further details concerning 
H.E.S.S. can be found in Hinton (2004) and references therein. 

H.E.S.S. observed Coma during the 2006 season for a total 
of 8.2 hr (corrected for the detector deadtime) comprising 19 
runs of duration ~28 min each. Those runs were accepted for 
data analysis if they met the quality control criteria described in 
Aharonian et al. !2004! Data were analysed using the moment- 
based Hillas analysis procedure described in Aharonian et al. 
(2006b ). Minimum cuts on the Cherenkov image sizf] of 80 
and 200 photoelectrons corresponding to standard and hard cuts 
were employed. The average zenith angle of the dataset was 
~53° yielding energy thresholds (peak detection rate for a power 
law source spectrum with an exponent of 2.1) of ~ 1.1 TeV and 
~2.3 TeV for standard and hard cuts analyses. This analysis 
follows on from preliminary H.E.S.S. results (Domainko et al. 
120071) . 

The large FoV of H.E.S.S. is well-suited to Coma as TeV 
emission could be expected from a variety of sites — the 
central radio halo or core; the radio-relic and adjacent galaxy 
merger/infall region; the degree-scale accretion shock suspected 
to surround the cluster (e.g. Voit 2005), and individual mem- 
ber galaxies. TeV y-ray significance skymaps covering a 7° x 7° 
FoV (from a mosaic of pointings) are presented in Fig.[T]), em- 
ploying oversampling radii of 0.2°, appropriate for moderately 
extended sources in the Coma field. Skymaps employing a 0.1° 
oversampling radius (Fig. !A.l! ) for pointlike sources are avail- 
able in the appendix. The CR background estimate in skymaps 
shown here is based on the template-model (Rowell 120031 ). em- 
ploying a region spatially overlapping the source region but not 
containing any y-ray-like events. Also available in the appendix 
are the distributions of skymap significances (Fig. !A.2b which 
are well-explained by Gaussians with standard deviation within 
a few percent of unity and means very close to zero, indicat- 
ing that the background estimate performs well over the FoV. 
Similar results were also obtained when employing alternative 
CR background estimates such as the ring/ring-segment and re- 
flected region models (Berge et al. 120071 ) which were used for 
upper limit calculations (in Table [!])• Results were also cross- 
checked using an alternative analysis chain. 

Table [T] summarises results for various locations in the 
Coma field guided by results from the ROSAT all sky survey 
(Briel et al. [1992] Voges et al. |1999] ), XMM-Newton observa- 
tions (Feretti & Neumann 120071 ) and Arecibo-DRAO radio ob- 
servations (Kronberg et al. 120071 ). The H.E.S.S. TeV excess sig- 

' Total photoelectron signal in the Cherenkov image. 
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Fig. 1. Left: Skymap of H.E.S.S TeV excess significance (colour-scale over ±4cr) calculated using Li & Ma ( I1983I I over a 7°x7° FoV, 
employing the template CR background model (Rowell'2003 ). An oversampling radius of 0.2° was used, appropriate for extended 
source searching. Overlaid contours (light-grey solid lines) represent total band (0.1 to 2.4 keV) smoothed X-ray counts s"' in 
log-scale from the ROSAT all-sky survey (Voges et al. |1999l l. Right: As for Left but with overlaid contours from radio observations 
(0.4 GHz - K contours rebinned from the original above 2.9K or ~11(t) from Kronberg et al. (|2007jl with strong point sources 
removed. The white dashed circle indicates the intrinsic 0.4° radius source size and position for the Coma Core (Tab.[T]) 



nificance S and flux upper limits O'^* (for an E"^ ' spectrum 
and using the method of Feldman & Cousins 09981 1) for E > I, 
> 5, and > 10 TeV were taken from standard, hard, and hard 
cuts analyses respectively. CR background estimates were taken 
from the reflected model (Berge et al. (i20(T7] l). 



In X-rays, extended emission from the Coma cluster and 
emission further to the southwest are evident. The southwest 
thermal X-ray emission is not entirely spatially coincident with 
the radio-relic (discussed below), but is centred on the galaxy 
sub-group NGC 4839 (labeled NGC 4839 in Table [B, ~ 20' 
closer to the Coma cluster core compared to the radio-relic. This 
sub-group is thought to represent a merger or infall of galaxies 
associated with Coma. Hard X-ray (18-30 keV) observations 
with INTEGRAL (Eckert et al. |2007 ) suggest excess emission 
in the direction of this infalling sub-group close to the cluster 
centre. In radio, the Coma core is visible as well as a radio ex- 
tension to the southwest known as the radio relic region (labeled 
1253-i-275-Relic in Table [1]). Two additional regions were also 
chosen to overlap the diffuse radio features labeled 'A' and 'B' 
from Kronberg et al. (I2007I I (at 0.4 GHz) for which circular re- 
gions of radii 0.9° and 0.5° respectively were used. Radio 'A' 
appears to well encompass the radio core of Coma which is dis- 
cussed at length in Thierbach et al. (2003) from their > 2GHz 
observations, whilst radio 'B' is a new feature from Kronberg 
et al. (I2007I I. In all cases no evidence for TeV emission was 
seen and 99% confidence level flux upper limits (assuming an 
E^^-^ spectrum) at several energy thresholds (E > 1, > 5 and 
> 10 TeV) were set. We note the highest excess significance fea- 
ture at ~ 4.1cr towards RA=12''55'" Dec=-i-27°15' is expected 
by chance given the number of independent trials (~ 10^) in the 
image. 



3. Discussion 

One of the most important properties of clusters of galaxies 
is the fact that CR protons remain diffusively confined in the 
magnetised intracluster medium for cosmological time scales. 
The maximum energy that can be confined depends on the (un- 
known) diffusion coefficient but an often made assumption is 
that the maximum energy is well above that relevant for TeV 
y-ray emission (Volk et al. 1 19961 . CR protons lose their en- 
ergy mainly via proton-proton interactions in the intergalactic 
medium. Due to the low density of this medium, the energy 
loss time is longer than the Hubble time. This implies that the 
hadronic CR content of a cluster is simply the superposition of 
the contributions from all the CR sources which have been ac- 
tive during the cluster lifetime, with little attenuation due to en- 
ergy losses. Under reasonable assumptions on the CR accelera- 
tion efficiency, the total non-thermal energy stored in the intra- 
cluster medium might be of the same order of magnitude of the 
thermal energy. For example Ryu et al. (2003) have estimated 
a non-thermal energy fraction reaching 50% of the thermal en- 
ergy. Such an amount of CR protons would result in copious 
emission of y-rays from the decay of neutral pions produced in 
proton-proton interactions. Since the most optimistic theoretical 
predictions are well within the capabilities of current-generation 
Cherenkov telescopes, the upper limits obtained by H.E.S.S. can 
be used effectively to constrain the non-thermal energy content 
of the Coma cluster. 

A remarkable feature of the y-ray emission from neutral pion 
decay is that its spatial profile is expected to follow the density 
profile of the gas which constitutes the target for proton-proton 
interactions. In the case of the Coma cluster, this gas is concen- 
trated within a core of radius ~ 300 kpc which, at the distance of 
the Coma cluster, coiTesponds to ~ 0.2° (see discussion below). 
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Table 1. Numerical summary for various regions in the Coma galaxy cluster and surrounding field. 



Name 


IRXS 


R.A. 
[J2000.0] 


Dec 
[J2000.0] 


'Rol rdeel 


-T [hi 


[0-] 

{E > 1,5, 10 TeV) 


Flux U.L. "♦O"'* 
{E >l,5,10TeV) 


Coma Core 


J 125947.9+275636 


12'' 59'" 


47.9' 


27°56'36" 


0.0 


7.3 


-0.5 


-1.2 


-1.4 


6.1 


0.3 


0.1 












0.2 


7.3 


+0.4 


-0.4 


-0.6 


10.8 


0.9 


0.5 












0.4 


7.3 


+ 1.1 


-0.5 


-1.5 


25.5 


1.7 


0.6 


1253+275-Relic§ 




12''55"' 


15.0' 


27°15'00" 


§ 


5.7 


+ 1.3 


+2.0 


+ 1.0 


15.9 


2.6 


1.2 


Radio-A 




12''55"' 


00.0' 


28°00'00" 


0.9 


4.6 


+2.0 


-2.2 


-2.0 


78.7 


2.3 


1.4 


Radio-B 




13'' 00"' 


00.0' 


30°15'00" 


0.5 


1.6 


+2.5 


+0.7 


+0.0 


77.8 


7.0 


3.4 


NGC 4839 


J125710.8+272426 


12''57'" 


24.3' 


27° 29' 52" 


0.2 


6.8 


+0.4 


+ 1.7 


+ 1.2 


9.0 


1.9 


1.0 












0.4 


6.8 


-1.8 


+0.3 


+0.1 


6.7 


2.0 


1.3 



1. Source Region of Interest (Rol) intrinsic radius. Actual radii used are convolved with the analysis PSF. A zero value here refers to a point-source analysis. 

2. Obsei'vation time (hr) coirected to a 0.7° off-axis angle using a standard cuts E > \ TeV response curve. 

3. Statistical significance using Li & Ma fT983^ 

4. 99% C.L. flux upper limit XIO"'^ ph cm"^ s"' 

§ Elliptical legion (0.33° X 0.2° with position angle 45°) as defined in Feretti & Neumann 420071 



This is the basis for the angular regions from which upper limits 
have been extracted. 

Table|2]demonstrates how the upper limits on the y-ray emis- 
sion convert into upper limits on the ratio of the cluster ther- 
mal energy to that of CR protons (non-thermal energy). This 
non-thermal to thermal energy ratio is denoted rf - EcRjEth- 
The Coma cluster thermal energy has been evaluated using the 
gas density profile and the intracluster medium temperature de- 
rived from X-ray data (e.g. Neumann et al. 2003) and resulted in 
E,h ~ 3.9 X lO*"^ and 1 .4 x lO*"^ erg for regions within 0.2 and 0.4 
degrees from the cluster centre respectively. The expected y-ray 
emission has been computed following Kelner et al. ( 120061) and 
assuming that the energy spectrum for CR protons is a single 
power law with spectral index a = 2.1 and 2.3 starting at an en- 
ergy of 1 GeV. The assumption of such hard spectra is justified 
by the fact that, due to CR confinement within the intracluster 
medium, the equilibrium spectrum must be equal to the CR in- 
jection spectrum at the sources. Note that much steeper spectra 
(a up to 6) are indicated for CRs accelerated at weak merger 
shocks (Gabici & Blasi 2003, Berrington & Dermer 2003"). The 
cluster non-thermal energy has been obtained by integrating the 
spectrum above 1 GeV and the resulting y-ray emission cor- 
rected for absorption in the cosmic infrared background (CIB) 
using the Salpeter initial mass function opacity given in Primack 
(1200 11 1. More recent constraints on the CIB from the TeV blazar 
lES 0229+200 (Aharonian et al. 2007 ) provide only a neghgible 
change in absorption from 1 to 10 TeV given Coma's proximity. 
Our upper limits can then be used to constrain this overall CR 
spectrum and hence non-thermal energy. 

The only missing piece of information is the spatial distri- 
bution of CRs. This quantity is not unambiguously known and 
it depends on the spatial distribution of CR sources in clusters. 
Here, two distinct situations were considered as two extreme 
cases. In the first one, referred to as Model-A, the radial profile 
of the CR energy density was assumed to follow the thermal en- 
ergy profile. In Model-B, a spatially homogeneous distribution 
for the CRs was assumed. 

In all cases considered, the most constraining data points in 
terms of energy flux are the upper limits for photon energies 
above 5 TeV. Our >10 TeV limits are marginally higher (~20%) 
whilst the > 1 TeV limits are factor 2 to 3 higher. For a CR spec- 
trum with a = 2.1 and assuming that the CR energy density 
follows the thermal energy density (Model-A) values for rf < 0.2 
for both the considered regions (0.2° and 0.4°) were obtained. 
In order to check how the upper limits depend on the assump- 



Table 2. Constraints ' on the ratio of CR (non-thermal) to thermal 
energy {EcR/E,f,) for the Coma cluster core region (within two 
radii) and assumed cosmic-ray distribution models A and B (see 
text). 



Radius 


a 


Model 


1 - Ecr/E,i, 


EcR [erg] 


0.2° (0.33 Mpc) 


2.1 


A 


<0.19 


<7.4xl0*' 


0.4° (0.67 Mpc) 


2.1 


A 


<0.18 


<2.5xl0" 


0.4° (0.67 Mpc) 


2.1 


B 


<0.25 


<3.5xlO*^2 


0.4° (0.67 Mpc) 


2.3 


A 


<0.55 


<7.7xl0« 



t The upper limit for E > 5 TeV has been used here. 



tions made on the CR spectrum and spatial distribution, the size 
of the region was fixed to 0.4° and two more cases were con- 
sidered: a homogeneous distribution of CRs (Model-B) and a 
softer CR spectrum with a - 2.3. The upper limits on rj are in 
these cases less stringent and are reported in the third and fourth 
rows of Table 12] resulting in < 0.25 and 0.55 respectively. In 
combination with our limits for £ > I TeV and >10 TeV, model- 
independent constraints on the ^ 10, ^ 50, ^100 TeV CR proton 
population in Coma were set. 

Our upper limits can be compared with those obtained from 
observations at other wavelengths and also other models. Firstly, 
our limits are slightly more constraining compared to those ob- 
tained from EGRET MeV/GeV data (?7=0.45 and 0.25 assuming 
a-2.\, 2.3; Pfrommer & EnBlin 2004). Our limits also rule out 
some of the models for CR acceleration in clusters of galaxies 
that predict high rj values, even up to 50% (Ryu et aI. I2003l l. The 
high frequency non-thermal radio emission of the Coma cluster 
(Thierbach et al. 120031 ) has been used by Reimer et al. ( 120041 ) to 
constrain rj by noticing that the radio emission from secondary 
electrons produced by CRs in proton-proton interactions cannot 
exceed the measured value. Recent observations (Brunetti et al. 
2008 ) of a steep radio spectrum from the cluster Abell 521 have 
also been considered in a similar way as in the Reimer et al. 
work on Coma. Their rather stringent limits obtained using ra- 
dio data {rf-lQ^^ to ~0.3) depend quadratically on the value of 
the intracluster Z?-field which has a large uncertainty of a fac- 
tor -10 {B =0.1 to 2juG). Additionally the < 10 GHz radio 
measurement constrains formally the < 0.1 TeV CR popula- 
tion within the range of B-fields used. ULs from H.E.S.S. and 
other VHE gamma-ray instruments make direct constraints on 
the £■ ^ 10 TeV CR population, energies well above that im- 
plied by Reimer et al. (|2004| |. and are essentially independent 
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of the B-field. Note that an additional preliminary E > 0.3 TeV 
upper limit from -19 hr of VERITAS (Perkins et al. 2008 ) obser- 
vations for the Coma core (0.3° radius) has been reported at ~3% 
Crab (or ~ 2 X lO''^ ph cm ^ s '). This VERITAS limit, when 
converted to an energy flux, provides a constraint on rj very sim- 
ilar to ours, albeit for CRs of slightly lower energies ^few TeV. 
Finally, our rj constraints for Coma are within a factor of two to 
three larger than those obtained from somewhat deeper H.E.S.S. 
observations (-33 hr) on the Abefl 85 cluster (M = 7.6 X 10^"^ 
Mq; z =0.055) (Ahai-onian et al. 2008). In this work, the con- 
straints 77 < 0.06 to 77 < 0.15 for Abell 85 were reported for a 
range of mass profiles and a CR spectral index of -2.1. 



4. Conclusions 

H.E.S.S. observed the Coma galaxy cluster for ~8 hr, obtaining 
upper limits to the E > 1, > 5, and > 10 TeV y-ray flux from the 
cluster core (0.2° and 0.4° radii). Additional regions for which 
upper limits are given include the radio-relic (1253-1-275); the 
merger/infall region associated with NGC 4839; and two large- 
scale radio features 'radio-A' (0.9° radius) and 'radio-B' (0.5° 
radius). 

Our results were compared to a model for the proton-proton 
y-ray emission assuming a proton spatial profile matching the 
centrally peaked thermal gas and injection spectral index of 2.1. 
In this case our £ > 5 TeV H.E.S.S. upper limit for the Coma 
core region within a 0.2° radius region (amounting to ~8% Crab 
flux units or ~ 10^'^ph cm"^ s"') constrains the fraction of en- 
ergy 77 in CRs to <0.2 times the thermal intracluster medium 
energy, or ~10''^ erg. This can be compared with the generally 
more stringent constraints on rj so far from radio observations 
(Reimer et al. 2004) in the range 77= lO"** to ~0.3 times the ther- 
mal energy for Z?-fields 0.1 to IjiG in the intracluster medium 
and proton injection spectral index a-2.\ to 2.5. It should be 
noted that the H.E.S.S. E > 5 TeV upper limit assuming a -2.1 
formally constrains ^ 50 TeV CRs whilst the radio limits per- 
tain to < 0.1 TeV CRs, highlighting the complementarity of 
the two approaches. The H.E.S.S. constraints for the amount 
of CRs stored in the intracluster medium also rule out the most 
optimistic theoretical models for CR acceleration in clusters of 
galaxies. 

Our upper limit for the NGC 4839 merger/infall group may 
also be useful, in conjunction with the hard X-ray emission, in 
constraining models for the additional shock acceleration of par- 
ticles potentially associated with this region. 

Our focus here has been on the expected centrally peaked 
proton-proton y-ray emission, but emission expected from the 
inverse-Compton scattering of electrons, and UHE proton/y in- 
teractions py — » e^e" + p' may follow the spatial profile of 
the degree-scale cluster accretion shock as an annulus of radius 
~ 1 - 2°. Inhomogeneities in the shock structure leading to TeV 
hotspots of up to degree in size (Keshet et al. 120031) . and uncer- 
tainties in the shock size prevented a specific attempt to search 
for such emission. TeV y-ray observations can in principle pro- 
vide direct constraints on the ability of such shocks to acceler- 
ate particles to multi-TeV energies and beyond. The significance 
skymap presented in Fig.[T]for 0.2° radii sources would suggest 
however that no such evidence for moderate-scale emission in 
the Coma field is present. Indeed, similar skymaps are obtained 
for a 0.4° radius. 

Overall our results indicate that deeper observations (to- 
wards 100 hr) of the Coma cluster, or other similarly massive and 
nearby galaxy cluster in TeV y-rays, are warranted to probe the 



Universe's largest-scale shocks. Such a detection appears pos- 
sible unless the total energy in the form of multi-TeV CRs is 
significantly less than ~ 10 % of the cluster thermal energy. 

Finally, the recently launched LAT instrument onboard the 
Fermi GST will provide critical constraints in the MeV/GeV 
band (likely within its first year or so of observation) and these 
results are eagerly awaited. 
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Appendix A: Additional Figures 

P.O. Box 103980, D 69029 H.E.S.S. images as for FiglUbut using a 0.1° integration radius 

are presented here as well as distributions of significances for 
both integration radii. 
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Fig. A.l. H.E.S.S. signficance skymaps as for Fig[T]but using an 
oversampling radius of 0.1°. Top: Overlaid contours (light-grey 
solid lines) represent total band (0. 1 to 2.4 keV) smoothed X-ray 
counts s"' in log-scale from the ROSAT all-sky survey (Voges 
et al. |1999l l. Bottom: Overlaid contours from radio observations 
from Kronberg et al. (120071) with strong point sources removed. 
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Fig.A.2. Distributions of H.E.S.S. significance from Figs: [T| 
and lA.ll for various oversampling radii — Top: 0.2°; Bottom: 
0.1°. The dashed line represents a Gaussian of zero mean and 
unit standard deviation. 



List of Objects 

'AGO 1656' on page|l 



